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Abstract 

Cosmic ray data may allow the determination of the proton-air cross 
section at ultra-high energy. For example, the distribution of the first 
interaction point in air showers reflects the particle production cross 
section. As it is not possible to observe the point of the first interaction 
Xi of a cosmic ray primary particle directly, other air shower observ- 
ables must be linked to Xi . This introduces an inherent dependence of 
the derived cross section on the general understanding and modeling 
of air showers and, therfore, on the hadronic interaction model used 
for the Monte Carlo simulation. We quantify the uncertainties arising 
from the model dependence by varying some characteristic features of 
high-energy hadron production. 

1 Introduction 

The natural beam of cosmic ray particles extends to energies far beyond the reach of any earth- 
based accelerator. Therefore cosmic ray data provides an unique opportunity to study interactions 
at extreme energies. Unfortunately, the cosmic ray flux is extremely small making direct mea- 
surements of the particles and their interactions impossible above ~ 100 TeV. One is forced to 
rely on indirect measurements such as extensive air shower studies, where interpretation of the 
data is very difficult. 

In this contribution we will briefly discuss different methods of measuring the proton-air 
cross section, focusing on methods that are based on extensive air shower (FAS) data. Figure [T] 
shows a compilation of proton-air cross section measurements and predictions of hadronic inter- 
action models currently used in cosmic ray studies [1-11] 

2 Methods of cross section measurements using cosmic ray data 
2.1 Primary cosmic ray proton flux 

Already in the 60's first estimates of the proton-air cross section cjp-air were made using cosmic 
ray data [1]. These early measurements are relying on two independent observations of the flux 
of primary cosmic ray protons after different amounts of traversed atmospheric matter. Firstly 
the primary proton flux <^{Xtop) is measured at the top of the atmosphere with a satellite or at 
least very high up in the atmosphere on a balloon at X^op = — 5 gcm^^. The second flux 
$(Xbottom) is measured with a ground based calorimeter at Xbottom = 600 — 1000 gcm~^, 
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Fig. 1: Current data of proton-air production cross section measurements [1-1 1] and model predictions [12-18]. 



preferentially at high altitude and using efficient veto detectors to select unaccompanied hadrons. 
The effective attenuation length can then be calculated straightforwardly from 

Aprod = (-^bottom — -'^top)/ log('I>top/*J'bottom)- (1) 

As it is impossible to veto aU hadronic interactions along the cosmic ray passage through the 
atmosphere, this attenuation length can only be used to obtain a lower bound to the high energy 
particle production cross section 

(m) 

'7p— air ^ T i (2) 
^prod 

where (m) is the mean mass of air. The method is limited to proton energies lower than '^TeV, 
since no sufficiently precise satellite or balloon borne data is available above this energy. By 
design the unaccompanied hadron flux is only sensitive to the particle production cross section, 
since primary protons with interactions without particle production cannot be separated from 
protons without any interaction. 



2.2 Extensive air showers 

In order to measure crp_air at even higher energies it is necessary to rely on EAS data [6-1 1]. The 
characteristics of the first few extremely high energy hadronic interactions during the startup of 
an EAS are paramount for the resulting air shower. Therefore it should be possible to relate EAS 
observations like the shower maximum Xmax. or the total number of electrons Nc{X) \x=x ^ ~ 
TVf^ and muons Nf^(X)\^_-^ — NJ^^'^ at a certain observation depth ^obs> to the depth of the 




Fig. 2: Definition of variables to characterize HAS longitudinal profiles. 



first interaction point and tlie cliaracteristics of the liigli energy liadronic interactions. 
Ground based observations 

In case of ground based extensive air shower arrays, the frequency of observing EAS of the 
same energy at a given stage of their development is used for the cross section measurement. By 
selecting EAS of the same energy but different directions, the point of the first interaction has to 
vary with the angle to observe the EAS at the same development stage. The selection of showers 
of constant energy and stage depends on the particular detector setup, but the typical requirement 
is [N^^^, NJ^'^'^^ = const at observation level. 

With the naming conventions given in Fig. |2l the probability of observing a shower of a 
given energy Eq and shower stage at the zenith angle 9 can be written as 



1 dN 



N dcose 



/r g-^l/->*int 
dAXi / dAX2 — 



X Pi(AXi) X P2{AX2) 

Here Xgtage defines the distance between the first interaction point and the depth at which the 
shower reaches a given number of muons and electrons as defined by the selection criteria. The 
experimentally inferred shower stage at observation level ^stage does, in general, not coincide 
with the true stage due to the limited detector and shower reconstruction resolution. This effect 
is accounted for by the factor Pres- The functions Pi and P2 describe the shower-to-shower 
fluctuations. The probability of a shower having its maximum at Xmax = -'^^i + AXi is expressed 
by Pi. The probability P2 is defined correspondingly with Xgtagc = AXi + AX2. 

In cross section analyses, Eq. © is approximated by an exponential function of sec 9. 
Assuming that the integration of Q over the distributions Pi, P2, and Pres does not yield any 
generally non-exponential tail at large sec 9, it can be written as 
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Fig. 3: Model predictions for the secondary particle multiplicity in high energy hadronic interactions. The lines denote 
the mean while the error bars indicate the RMS of the distributions. The shaded area is the range for SIBYLL using 

0-3 > /lOEeV > 3. 

However, the slope parameter A^^^ does not coincide with the interaction length Ajnt due to non- 
Gaussian fluctuations and a possible angle-dependent experimental resolution. Therefore the 
measured attenuation length can be written as 

S S 

^obs = ^int • kAXi ' kAX2 ' ^resolution = ^int • ^S- (5) 

The fe-factors k^Xi , ^AX2 ^rid A;^esoiution parametrize the contributions to A^^^ from the corre- 
sponding integrations. However, these integrations are difficult to perform separately and the 
individual /c-factors are not known in most analyses (for a partial exception, see [10]). 

Observations of the shower maximum X^ax 

Observing the position of the shower maximum directly allows one to simplify ^ by removing 
the term due to the shower development after the shower maximum P2. Also the detector reso- 
lution Pj-cs is much better under control for X^ax and can be well approximated by a Gaussian 
distribution. The resulting distribution is 

dXi / dAXi ^ ^ X Pi(AXi) X P,es(^max - ^max), (6) 

with Xi + AXi = Xmax- In analogy to Eq. (01) only the tail of P{X^^l^) at large X^^x is 
approximated by an exponential distribution 

whereas the exponential slope Aobs can be deduced from the convolution integral ([6]) as 

^^bs = '^int ■ ^AXi ' ^resolution = '^int ' ^X- (8) 
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Fig. 4: Mean and RMS values for the resulting Xmax, A''M(-'f)lx=iooo gcm-2 ^c(-'i^)|jf^iooogcm-2 distributions 
as a function of /ioecV using SIBYLL 2.1. For each data point with changed multiplicity, 1000 air showers are 
simulated and 10000 for a changed cross section. The dashed lines are polynomial fits of 2"^* order to guide the eye. 



Again /caxi and ^resolution contributions to A^^^ from the corresponding integrations 

of®. 

It was also recognized that Q can be unfolded directly to retrieve the original Xi -distribution, if 
the AXi -distribution is previously inferred by Monte-Carlo simulations [11]. Recently this trig- 
gered some discussion about the general shape and model dependence of the AXi -distribution 
[19]. This directly implies a corresponding model dependence of the A;aXi -factors. 



3 Impact of high energy interaction model characteristics on air shower development 

To explore the impact of uncertainties of the present high energy hadronic interaction models on 
the interpretation of EAS observables, we modified the CONEX [20] program to change some 
of the interaction characteristics during EAS simulation. To achieve this, individual hadronic 
interaction characteristics are altered by the energy-dependent factor 

f 1 E<1 PeV 

^^^^ "11 + (/lOEeV - 1) • logio(^/lPeV)/logio(10EeV/lPeV) E>1 PeV 

which was chosen to be 1 below 1 PeV, because at these energies accelerator data is available 
(Tevatron corresponds to 1.8 PeV). Above 1 PeV, f{E) increases logarithmically with energy, 
reaching the value of /ioeoV at 10 EeV 

The factor f{E) is then used to re-scale specific characteristic properties of the high energy 
hadronic interactions such as the interaction cross section, secondary particle multiplicity or in- 
elasticity. Obviously by doing this we may leave the parameter space allowed by the original 
model, but nevertheless one can get a clear impression of how the resulting EAS properties are 



depending on the specific interaction ciiaracteristics. 

We demonstrate tiie impact of a changing multiplicity nmuit and cross section a on the follow- 
ing, important air shower observables: shower maximum Xmax> and the total number of electrons 
N^'^'^, as well as muons NJ^'^'^ arriving at an observation level of X^hs = 1000 gcm~^. Figure [3] 
shows the range of extrapolations of rimuit used by the current hadronic interaction models and 
thus motivates the energy dependent re-scaling of nmuit by 0.3 > /loEeV ^ 3. 
All simulations are performed for primary protons at 10 EeV using the SIBYLL 2. 1 [17] interac- 
tion model. Figure |4] summarizes the results, which are discussed below. 

Multiplicity of secondary particle production 

The effect of a changed multiplicity on the Xmax-distribution is a shift to shallower Xmax with 
increasing rimuit- This is what is already predicted by the extended Heitler model [21] 

E 

Xmax oc A^, • In , (10) 

"'mult • -t/crit 

where is the electromagnetic radiation length and i^cri™ critical energy in air. This is 
a consequence of the distribution of the same energy onto a growing number of particles. The 
resulting lower energy electromagnetic sub-showers reach their maximum earlier. The impact 
on the RMS of the Xmax-distribution is small, but there is a trend to smaller fluctuations for an 
increasing number of secondaries. 

The total muon number after 1000 gcm~^ of shower development is rising if the multiplicity 
increases. This reflects the overall increased number of particles. The fluctuations are not signif- 
icantly affected. 

More interesting is the impact on the electron number A'^^'^'^, which shows a minimum close to 
/lOEeV = 1- The rising trend in the direction of smaller nmuit can be explained by the increase 
of Xjnax and therefore the shower maximum coming closer to the observation level. On the other 
hand the rising trend in the direction of larger nmuit is again just the consequence of a generally 
growing number of particles. In contr ary to the muon number the RMS does significantly change 
while rimuit gets larger. This can be explained by the strong dependence of fluctuations in A'^^'^*^ 
from the distance to the shower maximum. 

Cross section 

By construction, scaling the cross section does affect all hadronic interactions above 1 PeV, not 
only the first interaction. 

The mean as well as the RMS of the Xmax-distribution are decreasing with an increasing cross 
section. The effect is very pronounced, since the depth of the first interaction Xi is affected as 
well as the shower startup phase. Both effects are pointing to the same direction. This makes 
-'^max a very sensitive observable for a cross section measurement. 

The impact on the muon number NJ^^"^ is not very large. Since the shower maximum moves 
away from the observation level with increasing cross section, we just see the slow decrease of 
the muon number at late shower development stages, while the fluctuation of NJ^^^ stay basically 
constant. 

The mean electron number as well as its fluctuations depend strongly on the distance of X^ax 
from the observation level. Combined with the influence of the modified cross section on X^ax 
this explains well the strong decrease of the mean N^'^'^ as well as the RMS with increasing cross 
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Fig. 5: Simulated fcx -factors (fcx = A^j^g/Aint) for SIBYLL at 10 EeV with modified cross section and multiplicity. 
A^[,g is obtained by an exponential fit to the tail of the resulting CONEX Xmax-distribution. An ideal detector is 
assumed, hence kx = kAXi ■ The polynomial fits of 2"^^ order are only plotted to guide the eye. 



section. At very small cross sections the shower maximum comes very close to the observation 
level, which can be observed as a flattening in the mean N^'^'^ and the decrease of the fluctuations 
in N^'^'^ against the trend of increasing fluctuations of the position of the shower maximum itself. 



4 Summary 

All methods of EAS-based cross section measurements are very similar and thus suffer from the 
same limitations. 

• The values of all A;-factors must be retrieved from massive Monte-Carlo simulations. All 
analysis attempts so far have only calculated the combined factor of ks, respectively kx- 

• /c-factors depend on the resolution of the experiment and can therefore not be transferred 
simply to other experiments. 

• fcx-factors are inherently different from fcs-factors and can therefore not be transferred 
from an Xmax-tail analysis to that of ground based frequency attenuation or vice versa. 

• It cannot be disentangled whether a measurement of Aobs can be attributed to Aint entirely 
or at least partly to changed fluctuations in AXi and/or AX2. 

• Generally the Pi and P2 distributions have a complex shape and therefore the integra- 
tions of ^ and ^ to yield the approximations and ^ are leading to non-exponential 
contributions. 

• Any non-exponential contribution creates a strong dependence of the fitted Aobs on the 
chosen fitting range [22]. A strong non-exponential contribution makes the A;-f actor anal- 
ysis unusable. 



• It can be shown that the Pi(AXi)-distributions is very sensitive to changes of the high 
energy hadronic interaction characteristics and thus P{AX) = f{a, n^uit, ...) is a func- 
tion of a , Timuit and other high energy model parameters. Consequently this also makes 
the A;-factors depending on the high energy interaction characteristics k = f{a, rimuit, ...), 
which certainly must be considered for any cross section analysis. 

In Fig. |5] we show how the here presented simulations can be used to quantify the uncertainty 
caused in the fc-factors due to the dependence on nmuit to about it ~ 0.1 for a variation of the 
multiplicity by a factor from 0.3 up to 3. It is clear that even without considering the multiplicity 
as a possible source of uncertainty the d-dependence of the A;-factors certainly needs to be taken 
into account. Otherwise a systematic shift will be introduced into the resulting cjp-air, since part 
of the observed signal in Aobs is wrongly assigned to Aint, while in fact it must be attributed to 
k{a, rimuit, • • • ) [23]. This has not been considered in any EAS-based cTp_air measurement so far. 
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